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ABSTRACT 

We consider the unusual evolutionary state of the secondary star in Cygnus X-2. 
Spectroscopic data give a low mass (M2 ~ 0.5 — 0.7M Q ) and yet a large radius (R2 — 
7R Q ) and high luminosity (L 2 — 15OL ). We show that this star closely resembles a 
remnant of early massive Case B evolution, during which the neutron star ejected most 
of the ~ 3M0 transferred from the donor (initial mass Mn ~ 3.6M0) on its thermal 
time-scale ~ 10 6 yr. As the system is far too wide to result from common-envelope 
evolution, this strongly supports the idea that a neutron star efficiently ejects the 
excess inflow during super-Eddington mass transfer. Cygnus X-2 is unusual in having 
had an initial mass ratio q\ = M-a/Mi in a narrow critical range near ~ 2.6. 
Smaller lead to long-period systems with the former donor near the Hayashi line, 
and larger q\ to pulsar binaries with shorter periods and relatively massive white dwarf 
companions. The latter naturally explain the surprisingly large companion masses in 
several millisecond pulsar binaries. Systems like Cygnus X-2 may thus be an important 
channel for forming pulsar binaries. 
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1 INTRODUCTION 

Cygnus X-2 is a persistent X-ray binary with a long or- 
bital period (P = 9.84 d, Cowley, Crampton & Hutchings 
1979). The observation of unambiguous Type I X-ray bursts 
(Smale, 1998) shows that the accreting component is a neu- 
tron star rather than a black hole. The precise spectroscopic 
information found by Casares, Charles & Kuulkers (1998), 
and the parameters which can be derived from it, is sum- 
marized in Table 1. The mass ratio q = M2/M1 — 0.34 
implies that mass transfer widens the system, and is there- 
fore probably driven by expansion of the secondary star. 
Normally in long-period low-mass X-ray binaries (LMXBs) 
this occurs because of the nuclear evolution of a subgiant 
secondary along the Hayashi line, with typical effective tem- 
peratures T e ff,2 — 3000 — 4000 K. However Casares et al.'s 
observations show that this cannot be the case for Cygnus 
X-2. The secondary is in the Hertzsprung gap (spectral type 
A9 III): use of Roche geometry and the Stefan-Boltzmann 
law gives L 2 ~ 15OL with T cff , 2 ^ 7330 K (see Table 
1). Moreover the mass ratio q ~ 0.34, and the assump- 
tion that the primary is a neutron star and thus obeys 
Mi SSi 2M©, implies that the secondary has a low mass 
(M a = qMi <, 0.68M Q ). In contrast, an isolated A9 III 
star would have a mass of about 4M© . More recently Orosz 



& Kuulkers (1998) have modelled the ellipsoidal variations 
of the secondary and thereby derived a model-dependent in- 
clination of i = 62.5° ± 4° which translates into component 
masses (Mi = 1.78 ± O.23)M and (M 2 = 0.60 ± 0.13)M Q . 

In this paper we consider explanations for the unusual 
nature of the secondary in Cygnus X-2. We find only one 
viable possibility, namely that this star is currently close to 
the end of early massive Case B mass transfer, and thus 
that the neutron star has somehow managed to reject most 
of the mass (~ 3M©) transferred to it in the past. In support 
of this idea, we show that this type of evolution naturally 
explains the surprisingly large companion masses in several 
millisecond pulsar binaries. 



2 MODELS FOR CYGNUS X-2 

In this Section we consider four possible explanations for the 
unusual nature of the secondary in Cygnus X-2. We shall find 
that three of them are untenable, and thus concentrate on 
the fourth possibility. 
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Table 1. Observed and derived properties of Cygnus X-2. 

spectroscopic period P 9 d .844(±3) 
radial-velocity amplitude K2 88.0 ± 1.4 km s 
rotational velocity V2 sin i 34.2 ± 2.5 km s — 1 

spectral type Sp(2) A9 ± 2 III 

mass ratio q = M 2 /Mi 0.34 ± 0.04 

effective temperature T c ff 2 (7330 ± 320) K 
primary mass Mi J> (1.43 ± O.1O)M 

secondary mass M 2 ^ (0.49 ± O.O7)M 

< (0.68 ± O.O8)M 
. 

secondary radius R2 8.93R 

,2/3 

secondary luminosity L2 



(%) 

2O7.7L (i|) ;i/iS (fey 



2.1 A normal star at the onset of Case B mass 
transfer? 

The simplest explanation is that the position of the sec- 
ondary in the HR diagram is just that of a normal star 
crossing the Hertzsprung gap. Since such a star no longer 
burns hydrogen in the core, this is a massive Case B mass 
transfer as defined by Kippenhahn & Weigert (1967), here- 
after KW. Provided that the initial mass ratio qi Ss 1 the 
binary always expands on mass transfer, which occurs on 
a thermal time-scale. Kolb (1998) investigated this type of 
evolution systematically and found that the secondary's po- 
sition on the HR diagram is always close to that of a single 
star of the same instantaneous mass. For Cygnus X-2 the 
A9 III spectral type would require a current secondary mass 
M 2 ~ 4M , and thus a primary mass Mi = M 2 /q ~ 12M . 
This is far above the maximum mass for a neutron star and 
would require the primary to be a black hole, in complete 
contradiction to the observation of Type I X-ray bursts from 
Cygnus X-2 (Smale, 1998). 

We conclude that the secondary of Cygnus X-2 cannot 
be a normal star. Accordingly we must consider explanations 
in which it is far bigger and more luminous than expected 
from its estimated mass ~ (0.49 — O.68)M . 



2.2 A stripped subgiant? 



The type of Case B evolution described in subsection (2.1) 
above is known as 'early', in that mass transfer starts when 
the donor's envelope is still largely radiative rather than 
convective (as it would become as the star approached 
the Hayashi line), and 'massive', meaning that the helium 
core has a large enough mass that upon core contraction 
it does not become degenerate but instead ignites central 
helium burning. The corresponding minimum core mass is 
about O.35M , corresponding to a total ZAMS mass of 
2 — 2.5M (depending on the assumed degree of convective 
overshooting during central hydrogen burning). For lower 
initial masses we have 'low-mass' Case B (Kippenhahn, Kohl 
& Weigert 1967) . Here the donor's helium core becomes de- 
generate and the envelope is convective. After a possible 
early rapid mass transfer phase in which the binary mass ra- 
tio q = M2/M1 is reduced to Sa 1, (Bhattacharya & van den 
Heuvel, 1991; Kalogera & Webbink, 1996) the donor reaches 
the Hayashi line and mass transfer is driven by its nuclear 
expansion under hydrogen shell burning. The star remains 
on the Hayashi line, increasing its radius and the binary pe- 



riod as its core mass grows. Webbink, Rappaport & Savonije 
(1983) describe this type of 'stripped giant /subgiant' evo- 
lution in detail, and indeed fit Cygnus X-2 in this way. 
However there is a clear discrepancy between the observed 
effective temperature and that required for a Hayashi-line 
donor, which should be about 4100 K in this case. Web- 
bink et al. (1983) appeal to X-ray heating by the primary 
to raise the temperature to the observed 7330 K, but re- 
mark that since the heating only operates on the side of 
the secondary facing the primary, one would expect a much 
larger orbital modulation of the optical flux than actually 
observed (AV^bs — 0.3 mag) unless the orbital inclination i 
is low. Simple estimates (see the Appendix) show that such a 
small modulation would require i < 13.4°. However the mass 
function and mass ratio for Cygnus X-2 can be combined to 
show that Mi sin 3 i = (1.25±O.O9)M , so such small inclina- 
tions would imply Mi > 1OOM , again clearly incompatible 
with the very strong observational evidence for a neutron 
star primary. A still stronger argument can be constructed 
on the basis that the spectral type of the secondary is not 
observed to vary during the orbital cycle. 



2.3 A helium white dwarf undergoing a hydrogen 
shell flash? 

It is known that newly-born helium white dwarfs can un- 
dergo one or more hydrogen shell flashes during their evo- 
lution from the giant branch to the white dwarf cooling 
sequence. During these flashes the star has a much larger 
photosphere. Calculations by Driebe et al. (1998) show that 
only low-mass He white dwarfs in the interval O.21M Ss 
Mwd Si O.3OM can undergo such a flash which in turn can 
put the star in the same position on the HR diagram as the 
secondary of Cygnus X-2 (e.g. the first shell flash of the se- 
quence with Mwd = O.259M ). However, the required low 
secondary mass has a price: the observed mass ratio implies 
Mi = qM 2 ~ (0.76 ± O.O9)M , much smaller than required 



by the mass function (Mi sin 



(1.25 ± O.O9)M ). To 



makes matters worse, the evolutionary track crosses the rel- 
evant region of the HR diagram in an extremely short time: 
the star's radius expands on a time-scale r = di/dlni?2 — 
33 yr. Not only does this give the present system an implau- 
sibly short lifetime, the radius expansion would drive mass 
transfer at a rate ~ AMh/t ~ few xlO _4 M yr -1 , again 
totally inconsistent with observations. We conclude that the 
secondary of Cygnus X-2 cannot be a low-mass He white 
dwarf undergoing a hydrogen shell flash. 



2.4 A star near the end of early massive Case B 
mass transfer? 



We saw in subsection (2.1) above that a secondary near the 



onset of early massive Case B mass transfer (i.e. with q Ss 1 
throughout) is ruled out for Cygnus X-2, as the required 
stellar masses conflict with observation. However, a more 
promising assignment is a secondary near the end of an early 
massive Case B evolution which began with qi JJ; 1. 

KW have investigated this process in detail. In contrast 
to t he case qi Ss 1 discussed by Kolb (1998), and considered 



in (2.1) above, the ratio q\ J; 1 means that the binary and 
Roche lobe initially shrink on mass transfer. Adiabatic sta- 
bility is nevertheless ensured because the secondary's deep 
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radiative envelope ('early' Case B) contracts on rapid mass 
loss. Mass transfer is therefore driven by the thermal-time- 
scale expansion of the envelope, but is more rapid than for 
gi ^ 1 because of the orbital shrinkage. Once M 2 is reduced 
to the point that q ^ 1, the Roche lobe begins to expand. 
This slows the mass transfer, and shuts it off entirely when 
the lobe reaches the thermal-equilibrium radius of the sec- 
ondary, since the latter then has no tendency to expand fur- 
ther (except possibly on a much longer nuclear time-scale). 
Calculations by KW and Giannone, Kohl & Weigert (1968), 
hereafter GKW, show that in some cases the orbit can shrink 
so much that the process ends in the complete exhaustion of 
the donor's hydrogen envelope, ultimately leaving the core 
of the secondary in a detached binary. Alternatively, the 
rapid Case B mass transfer may end with the donor on the 
Hayashi line, still retaining a large fraction of its original 
hydrogen envelope. However, there will be no long-lasting 
phase of mass transfer with the donor on the Hayashi line 
because the star starts shrinking with ignition of central he- 
lium burning (KW). Neither of these two cases describes 
Cygnus X-2. However, there is an intermediate possibility: 
the initial mass ratio gi may be such that the donor retains 
a small but non-negligible hydrogen envelope as mass trans- 
fer slows. The current effective temperature of 7330 K shows 
that the companion's envelope is mainly radiative, with only 
a very thin surface convection zone. (A paper in preparation 
by Kolb et al. shows this in detail.) In the example computed 
by KW the donor, at the end of mass transfer, is not on the 
Hayashi line, but at almost the same point in the HR dia- 
gram as it occupied immediately before mass transfer began 
(this is confirmed by the detailed numerical calculations of 
Kolb et al.). Just before the process ends we then have an 
expanding low-mass donor, driving a modest mass transfer 
rate in a long-period expanding binary, but at the HR dia- 
gram position of a much more massive normal star. As we 
shall show below, for an initial donor mass of about 3.6Mq 
the end point of such an evolution can be made to match 
closely that observed for the secondary of Cygnus X-2. (Note 
that we have not performed detailed numerical calculations 
for this paper, but rather used the results of KW and GKW. 
The forthcoming paper by Kolb et al. reports detailed cal- 
culations.) 

Clearly this idea offers a promising explanation of the 
secondary in Cygnus X-2. However there is an obvious dif- 
ficulty in accepting it immediately. KW's calculations as- 
sumed that the total binary mass and angular momentum 
were conserved, and in particular that the primary retained 
all the mass transferred to it. But the primary of Cygnus 
X-2 is known to be a neutron star, with a mass presumably 
2M0, so we must require instead that it accretes rela- 
tively little during mass transfer. This agrees with the idea 
that a neutron star cannot accrete at rates greatly in excess 
of its Eddington limit (~ 10~ 8 Mq yr _1 ), and the fact that 
almost all of the mass is transferred at much higher rates 
(k, 10" 6 M Q yr" 1 ). We thus follow earlier authors (Bhat- 
tacharya & van den Heuvel, 1991; Kalogera & Webbink, 
1996) in postulating that the neutron star is extremely effi- 
cient in ejecting most of the super-Eddington mass transfer, 
rather than allowing the excess mass to build up into a com- 
mon envelope. Clearly common- envelope evolution could not 
produce Cygnus X-2: the current binary period of 9.8 d 
shows that far too little orbital energy could have been re- 



leased to remove the envelope of any plausible progenitor 
for the secondary (see e.g. Section 4 below). By contrast, it 
is at least energetically possible to expel most of a super- 
Eddington mass transfer rate, provided that this is done 
at large enough distance R c j from the neutron star. If the 
matter is given just the escape velocity at _R e j the ratio of 
ejection to accretion rate is 



Ma, 



Rcj 

R, ' 



(1) 



where R* is the radius of the neutron star. Ejecting all but 
about 1% of the transferred matter thus requires 



R cj £ 10OR, ~ 10 8 cm 



(2) 



which is far smaller than the size of the accretion disc around 
the neutron star for example. (This point is discussed further 
by King and Begelman, 1999.) 

In the next section we will consider an early massive 
Case B evolution for Cygnus X-2. We will find that the hy- 
pothesis of efficient mass ejection by the neutron star allows 
excellent agreement with the current state of the system, 
as well as plausible explanations for the observed states of 
several detached pulsar binaries. 



3 EARLY MASSIVE CASE B EVOLUTION 
FOR NEUTRON-STAR BINARIES 

The main features of early massive Case B evolution can be 
understood by considering the relative expansion or contrac- 
tion of the donor's Roche lobe 7?l and the thermal equilib- 
rium radius _R 2jC which the donor attains at the end of mass 
transfer (see e.g. GKW). As discussed above, we assume that 
the neutron star ejects any super-Eddington mass inflow. 
Since the mass transfer rate exceeds the Eddington limit by 
factors J; 100 (see above), almost all of the transferred mass 
must be ejected, and to an excellent approximation we can 
assume that the neutron star mass Mi remains fixed dur- 
ing the mass transfer (the equation for J?l can actually be 
integrated exactly even without this assumption, but at the 
cost of some algebraic complexity) . We assume further that 
the ejected mass carries the specific angular momentum of 
the neutron star's orbit. This is very reasonable, since the 
ejection region is much smaller than the size of the disc (see 
eq. ^). Then we can use the result quoted by Kalogera & 
Webbink (1996) to write 



Rh,i 



M 2 



_ / M2A '' ( M\ e 2(M 2 -M m )/Mi 



4/3 



M ) 



(3) 



where M2 and M = Mi + M 2 are the donor and total binary 
mass at any instant, and M%i, Mi their values at the onset 
of mass transfer. In writing (^) we have used the simple 
approximation Rl/cl tx (M2/M) 1 / 3 , where a is the binary 
separation. Using this and Kepler's law we get the change 
of binary period P as 



P_ 
~P~i = 

so that 
P_ 



Rl 
Rla 



Mai 
M 2 



3/2 



M 2i 
M 2 



M 
~M 



1/2 



o 3(M 2 -M 2i )/M 1 



(4) 



(•») 
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Figure 1. Schematic mass radius diagram of early massive Case 
B evolution showing different outcomes of mass transfer depend- 
ing on the initial mass ratio and initial separation. Full line: 
thermal equilibrium radius i?2,c of a star with hydrogen shell 
burning and a non-degenerate He core of mass Mh c as a function 
of total mass M2 (generalized main sequence). i?HL is the maxi- 
mum radius attained at the Hayashi line, Rsc the radius at the 
Schonberg— Chandrasekhar limit. Dashed lines: Roche lobe radius 
i?L as a function of M2 for three different values of the initial mass 
ratio q\ and the maximum possible orbital separation. The dashed 
lines can be shifted downwards by an arbitrary amount subject 
to the condition that i?2i > ^SC- Dash-dotted line: mass trans- 
fer starting when the secondary has just reached the Schonberg- 
Chandrasckhar limit. This example shows that q cr ; t as defined in 
section 3 depends on the initial orbital separation of the binary. 



Conventional massive Case B evolution always begins with 
a mass ratio qi = M2i /Mi sufficiently large that i?L initially 
shrinks. From (^) it is easy to show that this requires qi > 
1.2. If during the evolution M2 decreases enough that q < 
1.2, J?l begins to expand again. The curves of logRh thus 
have the generic U-shaped forms shown in Figs. 1-3. 

The thermal equilibrium radius R2.C depends on the rel- 
ative mass M-He/M2 of the donor's helium core. Figs. 1-3 
show the so-called 'generalized main sequences' of GKW, the 
first schematically, and the latter two for M2i = 3M© , 5M0 . 
For a core-envelope structure to be applicable, the star must 
have at least finished central hydrogen burning. Since the 
mass transfer takes place on the thermal time-scale of the 
donor there is little nuclear evolution, and we can regard 
Mne as fixed. The quantity i?2,e shown in Figs. 1-3 there- 
fore gives the thermal equilibrium radius attained by the star 
after transferring varying amounts of its hydrogen envelope. 
The evolution of the system is now specified by the initial 
mass ratio qi and the radius i?2i of the donor at the onset of 
mass transfer. This can lie between the maximum radius Rb 
reached during central hydrogen burning and one almost as 
large as the value Rhl at the Hayashi line (mass transfer 
is adiabatically unstable if the donor develops a deep con- 
vective envelope). The allowed initial radius range is about 
a factor 2 for a donor with M2i = 2.5M0, increasing to a 
factor ~ 6 for M 2i = 5M© (Bressan et al. 1993). 

We see from the schematic Fig. 1 that three qualita- 
tively different outcomes of early massive Case B mass trans- 
fer are possible: 



1. 'small' qi, i.e. JVfe only slightly larger than Mi. Here 
the Roche lobe soon begins to expand, so the curves of log Rl 
and log R2.C cross before much mass is transferred. Mass 
transfer ends with ignition of central helium burning as this 
makes the star shrink somewhat. At this point the secondary 
still has a thick hydrogen envelope and lies on the Hayashi 
line, with the binary having a long orbital period. During 
central helium burning the secondary stays on the Hayashi 
line but remains detached. Mass transfer resumes only after 
central helium burning when the star approaches the asymp- 
totic giant branch (AGB). Mass transfer is again driven by 
nuclear evolution (double shell-burning). Since core helium- 
burning has ceased this is Case C. Here (unusually) the mass 
transfer is adiabatically stable, despite the secondary's deep 
convective envelope, since the mass ratio is already 1. The 
average mass transfer rate is again of order 10 _6 Mq yr~\ 
but may become up to an order of magnitude higher during 
thermal pulses (Pastetter & Ritter 1989). The donor may 
also lose large amounts of envelope mass in a wind. If the 
binary again manages to avoid common-envelope evolution 
by ejecting most of the transferred mass, mass transfer will 
finally end once the secondary's envelope has been lost, leav- 
ing a very wide (period ~ 10 yr) binary containing a neutron 
star and a CO white dwarf. 

2. 'critical' qi. We define this as the case where the 
Rl, R2,c curves cross at the 'knee' in the mass-radius curve, 
i.e. with M2 only slightly larger than Mh c , so the mass trans- 
fer depletes almost the entire envelope. The remnant donor 
retains a thin hydrogen envelope, and lies between the main 
sequence and the Hayashi line on the HR diagram. Thus the 
envelope mass is low enough to prevent the star lying on the 
Hayashi line, but not so low that the remnant is small and 
hot, i.e. to the left of the main sequence. The initial separa- 
tion must be small enough that mass transfer starts before 
central helium burning, but large enough that it starts only 
after the donor has reached the Schonberg-Chandrasekhar 
limit. This limit is defined as the point where the isothermal 
helium core has reached the maximum mass which is able 
to support the overlying layers of the star, i.e. the point at 
which core collapse begins. The orbital period is shorter than 
in case 1., but longer than in case 3. below. Cygnus X-2 is 
an example of this evolution, viewed at the point where the 
donor has almost attained its thermal-equilibrium radius, 
and mass transfer is well below the maximum thermal-time- 
scale rate. This evolution ends with nuclear evolution of the 
donor to smaller radii as the mass of the hydrogen-rich en- 
velope is further reduced by shell burning. The system de- 
taches, leaving a helium-star remnant which subsequently 
ignites central helium burning and finally becomes a CO 
white dwarf. 

3. 'large' qi. The curves cross only when the hydrogen 
envelope is effectively exhausted. The remnant is a helium 
star and the orbital period is short. If Mh c 5s O.9M0 the 
helium star evolves directly into a CO white dwarf. If Afee <; 
IMq this star re-expands during helium shell-burning. This 
in turn can give rise to a further phase of (so-called Case 
BB) mass transfer, e.g. Delgado & Thomas, 1981; Law & 
Ritter, 1983; Habets, 1985, 1986). 

Table 2 summarizes the possible Case B evolutions with 
a neutron-star primary. The various outcomes all reflect the 
general tendency of larger qi (i.e. larger M2i) to produce 
greater orbital contraction, and thus smaller and relatively 
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Table 2. Outcomes of Case B evolution with a neutron-star primary 

subcase intermediate stages Pf (d) final WD companion 

low mass Hayashi-line LMXB ~ 10 — 1000 He, obeys mass-period relation 

early massive, < (jr cr ; t Hertzsprung-gap XRB, Case C mass transfer <; 1000 CO, obeys mass-period relation 

early massive, qi ~ g cr ; t Cygnus X-2 ~ 10 CO, overmassive 

early massive, qi > q CI n NS + He star, Case BB mass transfer? Ss 1 — 10, CO, overmassive 




Mo/M, 



Figure 2. Mass-radius diagram for early massive Case B evo- 
lution. Full (and long dashed) line: generalized main sequence 
for 3Mq stars taken from GKW. Short dashed lines: Roche 
lobe radius R^ as a function of M 2 computed from (3) with 
Mi = 1.4M ,M 2 i = 3M for two values of R 2i (0.75fi H L and 
0.55R HL ). 



less massive remnants (their mass is a smaller fraction of a 
larger M 2 i) and short orbital periods. Conversely, larger R 2 i 
in the range Rb — Rhl, where Rb is the maximum radius 
reached during central hydrogen burning, produces exactly 
the opposite trends. To model Cygnus X-2 one thus needs 
qi close to the critical value. 

As can be seen from the following arguments the range 
of possible solutions is strongly constrained by equation (J3J) : 
with the assumption that Mi remains essentially constant 
during the evolution, M%i is fixed by the observed mass ratio 
qt = 0.34±0.04. On the other hand, qi is also essentially fixed 
by the model assumption that the current donor is close 
to the end of Case B mass transfer. This in turn means 
that the donor is now close to thermal equilibrium, with 
its luminosity therefore coming mainly from hydrogen shell 
burning. 

Use of the generalized main sequence for 3Mq stars 
given in GKW (their fig. 7, shown here in Fig. 2) demon- 
strates that one can account for the current state of the 
donor in Cygnus X-2 if this star had M 2 i — 3M (i.e. 
qi ~ 2.1) and 0.55-Rhl Si R21 0.75-Rhl when mass transfer 
began, corresponding roughly to 0.5Mq <S M 2f Si 0.7Mq. 
The numerical results of GKW (their Table 2 and Fig. 4) 
suggest that the current state of the donor is then given 
approximately by M 2 = O.79M ,M Hc = O.67M and thus 
go = M He /M 2 = 0.85; with T cff , 2 = 7060 K, L 2 = 126L 
and R2 = 7.5R . These values correspond to an orbital 
period P ~ 8.4 d. From (^) we infer that mass transfer 
started at an orbital period Pi ~ 3.5 d. These quantities are 




Figure 3. Mass-radius diagram for early massive Case B evo- 
lution. Full line: generalized main sequence for 5M stars taken 
from GKW. Short dashed line: Roche lobe radius R^ as a function 
of M 2 computed from (3) with Mi = 1.4M , M 2i = 5M 



close to those given in Table 1, although the predicted M 2 
is slightly larger than the observational estimate. This may 
either result from the fact that the chemical composition of 
the models defining the generalized main sequences is char- 
acterized by a step function at M r = M c and thus differs 
from that of evolutionary models, or from the fact that the 
initial chemical composition (X — 0.602, Z = 0.044) and the 
opacities used by GKW are rather outdated. In addition the 
core mass A^hc = O.67M inferred above points to an ini- 
tial mass higher than the ~ 3M suggested earlier. In fact 
to fit the observed value of O.5M ^ M 2 f ^ O.7M , more 
modern calculations than those of GKW (e.g. Bressan et al. 
1993) yield the required core mass if the initial mass was 
3.2 <, M 2i ^ 4.1M . 

We conclude that the observational data for Cygnus X- 
2 are well reproduced if we assume it is a remnant of early 
massive Case B evolution with qi close to the critical value 
~ 2.3 - 2.9, if Mi ~ 1.4M , or ~ 1.8 - 2.3, if we adopt the 
value Mi ~ 1.8M derived by Orosz & Kuulkers (1998). For 
the remainder of this paper we shall adopt AIi = 1.4M . 



4 END PRODUCTS 

The discussion above shows that for both low-mass Case B, 
and for early massive Case B with gi below a critical value 
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Table 3. Binary pulsars with relatively massive WD companions 



pulsar 


P(d) 


M 2 (M s ) 


P B (ms) 


J0621+1002 


8.32 


> 0.45 


28.8 


J1022+1001 


7.81 


> 0.73 


16.5 


J2 145-0750 


6.84 


> 0.43 


16.1 


B0655+64 


1.03 


0.7 : 


196 



(~ 2.6 for M 2 i ~ 3.6Mq), the evolution leads to a long- 
period binary with the donor on the Hayashi line. As is well 
known, the luminosity and radius of such a star are fixed 
by its degenerate core mass M c rather than its total mass 
M 2 . Even though the degenerate core is different in nature 
in the two cases, it is possible to give a single formula for 
the radius, i.e 

_ 3.7 x lQ 3 mj 

V2 l + m 3 + L75m 4- W 

where r 2 = R 2 /KQ,m c = M c /Mq (Joss, Rappaport & 
Lewis, 1987). Then using the well-known relation 

3/2 

P = 0.38^- d (7) 

m 2 ' 

(m 2 = M 2 /Mq) which follows from Roche geometry, we get 
a relation between M C ,M 2 and the orbital period P (e.g. 
King, 1988). Once all of the envelope mass has been trans- 
ferred we are left with a wide binary containing a millisecond 
pulsar (the spun-up neutron star) in a circular orbit with 
the white-dwarf core of the donor. Since at the end of mass 
transfer we obviously have m 2 = m c , such systems should 
obey the relation 

P~8.5xl0 4 f^ 3 m \* iwrt ) d - (8) 

The timing orbit of the millisecond pulsar allows constraints 
on the companion mass, so this relation can be tested by ob- 
servation. Lorimer et al. (1995), Rappaport et al. (1995) and 
Burderi, King & Wynn (1996) show that while the relation is 
consistent with the data for a majority of the ~ 25 relevant 
systems, there are several systems (currently 3 or 4; see Fig. 
4 and Table 3) for which the white dwarf mass is probably 
too large to fit. While one might possibly exclude B0655+64 
because of its long spin period P s (but see below) , the other 
three systems are clearly genuine millisecond pulsars. 

Our considerations here offer a simple explanation for 
this discrepancy. If the initial mass ratio qi lies above the 
critical value (~ 2.6 for M 2 i ~ 3.6M0), the donor radius 
will be less than Phl at the end of early massive Case 
B mass transfer, and the orbital period relatively shorter. 
When such systems finally detach from the Roche lobe, the 
WD companion is considerably more massive than expected 
for the orbital period on the basis of the Hayashi-line re- 
lation (|^). We see from (|E|) that systems with large initial 
companion masses M 2 [ 4M0) can end as short-period 
systems. Table 4 and Fig. 4 show the expected minimum 
final periods Pf and companion masses M 2 for various M 2 \, 
i.e. assuming that mass transfer began with the donor at the 
Schonberg-Chandrasekhar limit. 

Table 4 and Fig. 4 show that early massive Case B evo- 
lution with qi larger than the critical value can explain the 



Table 4. Short-period end-products of early massive Case B evo- 
lution 



M 2i (M Q ) 


Pf(min) (d) 


M 2 (M ) 


4.0 


6.60 


0.56 


5.0 


1.25 


0.79 


6.0 


0.268 


1.05 


7.0 


0.048 


1.46 



The minimum final period Pf(min) is calculated assuming that 
mass transfer began at the Schonberg-Chandrasekhar limit. Ini- 
tial donor masses M 2 \ 4Mq may possibly be ruled out because 
of delayed dynamical instability (see text). 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 

M WD /M G 



Figure 4. P — M 2 diagram for long-orbital-period binary ms- 
pulsars and Cygnus X-2. Full lines: predicted P — M 2 relation 
for systems in which the donor star evolved along the Hayashi 
line (taken from Rappaport et al. 1995). Dashed line: predicted 
minimum orbital period for systems in which the white dwarf 
secondary was formed via early massive Case B mass transfer. The 
short-dashed part of this curve is realised only if the neutron star 
can expel the transferred mass despite the binary being subject 
to the delayed dynamical mass transfer instability. Also shown 
are the (P, M 2 ) values of the four 'discrepant' ms-pulsar binaries 
listed in Table 3, and the current position of Cygnus X-2. 

discrepant pulsar systems of Table 3. Indeed it appears that 
the process can end with very short orbital periods, offer- 
ing an alternative to the usual assumption of a common- 
envelope phase (Bhattacharya, 1996; Tauris, 1996) The lim- 
iting factor for this kind of evolution may be the so-called 
'delayed dynamical instability' (Webbink, 1977; Hjellming, 
1989). For a sufficiently massive initial donor, mass transfer 
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eventually becomes dynamically unstable because the adi- 
abatic mass radius exponent of a strongly stripped radia- 
tive star becomes negative and the donor begins to expand 
adiabatically in response to mass loss. For typical neutron 
star masses Mi ~ 1.4Mq this would limit M 2 i to values 
<> (4-4.5)M (Hjellming 1989; Kalogera & Webbink, 1996). 
Given the uncertainties in our estimates, this is probably 
consistent with the initial mass M 2 \ J> 5M needed to ex- 
plain the current state of the most extreme discrepant sys- 
tem (B0655+64). The other three systems can all be fitted 
with M 2 \ Si 4M , so there is no necessary conflict with the 
mass limits for the delayed dynamical instability. 

We may now ask about the end states of such an evolu- 
tion if the neutron star were unable to eject the mass trans- 
ferred at the very high rates expected once the delayed dy- 
namical instability set in and the system went through a 
common envelope phase instead. Using the standard pre- 
scription for estimating the parameters of a post common 
envelope system (e.g. Webbink 1984) with the values for 
M 2i and M 2( = M 2 given in Table 4, Mi = 1.4M and 
QceA = 0.5, where A ~ 0.5 is a structural parameter and 
qce the common envelope efficiency parameter defined by 
Webbink (1984), we find that the final orbital period is 
0.04 ^ Pf(d) Si 0.2 if mass transfer set in when the donor 
was already near the Hayashi line, and smaller still if mass 
transfer set in earlier or if qce is smaller than unity. On the 
other hand, unless qce is significantly smaller than unity, 
common envelope evolution starting from a system with the 
donor on the asymptotic giant branch would end with peri- 
ods much longer than those of the systems listed in Table 3 
and shown in Fig. 4. Table 5 shows the outcome of common- 
envelope evolution in the two extreme cases where mass 
transfer starts at the Schonberg-Chandrasekhar limit, and 
on the Hayashi line. Although common-envelope efficiencies 
exceeding unity are discussed in the literature (i.e. energy 
sources other than the orbit are used to expel the envelope), 
the values found for /3ce = A qce in Table 5 show that 
to produce Cyg X-2-like systems requires absurdly large 
efficiencies, even starting from the most distended donor 
possible. Thus common envelope evolution does not offer 
a promising explanation for these systems. Their very exis- 
tence may thus indicate that an accreting neutron star can 
eject mass efficiently even at the very high mass transfer 
rates encountered in the delayed dynamical instability. We 
conclude therefore that even very rapid mass transfer on 
to a neutron star does not necessarily result in a common 
envelope (cf King & Begelman, 1999). 

We note finally that all of the pulsars of Table 3 have 
spin periods much longer than their likely equilibrium peri- 
ods (i.e. they lie far from the 'spinup line', cf Bhattacharya & 
van den Heuvel, 1991), suggesting that they have accreted 
very little mass (<< O.1M ) during their evolution. This 
agrees with our proposal that these systems are the direct 
outcome of a super-Eddington mass transfer phase in which 
almost all the transferred mass is ejected. 

5 SPACE VELOCITY AND POSITION OF 
CYGNUS X-2 IN THE GALAXY 

The distance to Cygnus X-2 derived from the observations 
of Type 1 X-ray bursts (Smale 1989) is d = (11.6 ±0.3) kpc. 



From the galactic coordinates I — 87.33° and b = —11.32° 
and the solar galactocentric distance R a = (8.7 ± 0.6) kpc 
one derives a galactocentric distance for Cygnus X-2 of 
dGC = (14.2 ± 0.4) kpc and a distance from the galactic 
plane of z = (-2.28 ± 0.06) kpc. Thus Cygnus X-2 has a 
very peculiar position indeed, being not only in the halo but 
also in the very outskirts of our galaxy. But not only is its po- 
sition peculiar, its space velocity with respect to the galactic 
centre is also surprising. It can be shown that the observed 
heliocentric radial velocity of 7 = (—208.6 ± 0.8) km s _1 
(Casares et al. 1998) is totally incompatible with prograde 
rotation on a circular, even inclined orbit around the galactic 
centre (Kolb et al., in preparation). The orbit is either highly 
eccentric and/or retrograde. In either case Cygnus X-2 must 
have undergone a major kick in the past, presumably when 
the neutron star was formed in a Type II supernova. Since 
prior to the supernova explosion the primary was much more 
massive (Mn ^ 1OM ) than the secondary (M 2i ~ 3.6M ), 
the latter was still on the main sequence when the super- 
nova exploded. Thus the age of Cygnus X-2 (and the time 
elapsed since the supernova) is well approximated by the 
nuclear time-scale of the secondary, which is ~ 4 10 8 yr for 
a ~ 3.6M star. This means that Cygnus X-2 must have 
gone around the galactic center a few times since its birth 
or supernova explosion and that, therefore, its birthplace in 
the galaxy cannot be inferred from its current position and 
velocity. 



6 DISCUSSION 

We have shown that the unusual nature of the secondary 
star in Cygnus X-2 can be understood if the system is near 
the end of a phase of early massive Case B evolution in which 
almost all of the transferred material is ejected. The system 
is unusual in having had an initial mass ratio qi = M2/M1 in 
a narrow critical range near qi ~ 2.6; smaller ratios lead to 
detached systems with the secondary near the Hayashi line, 
and larger ratios produce binary pulsars with fairly short 
orbital periods and relatively massive white dwarf compan- 
ions. During this evolution, much of the companion's origi- 
nal mass (~ 3M for Cygnus X-2) is transferred and conse- 
quently lost on the thermal time-scale ~ 10 6 yr of this star. 
Evidently the huge mass loss rate and the short duration 
of this phase make it difficult to detect any systems in this 
state; they would probably resemble Wolf-Rayet stars of the 
WNe type (i.e. showing hydrogen). 

Cygnus X-2 is currently near the end of the thermal 
time scale mass transfer phase, so that its mass transfer rate 
is now well below the thermal time-scale value, and proba- 
bly given by the accretion rate. At M acc ~ 2x 1O _8 M yr -1 
(Smale, 1998), this is nevertheless one of the highest in any 
LMXB, making it easily detectable. Only a full calculation 
of the evolution, with in particular a detailed model for the 
secondary, can predict the duration of the current phase; 
this is not any easy task, as this star deviates strongly from 
thermal equilibrium during most of the evolution. But it is 
clear that the mass transfer rate will decline as the remain- 
ing few tenths of a solar mass in the hydrogen envelope are 
transferred. Cygnus X-2's long orbital period and large ac- 
cretion disc mean that even its current mass transfer rate 
only slightly exceeds the critical value required for a persis- 
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Table 5. Outcomes of common-envelope evolution for early massive Case B 



Case 




M 2i (M s ) 


M 2{ (M Q ) 


Mi (M ) 


^2i (R©) 


P f (d) for /3 C e = 0.5 


/3 CE for P f = 8 d 


/3 CE for P f = 1 d 


Bo- — 


Ran 


4.0 


0.56 


1.4 


8.5 


0035 


86.9 


21.7 






5.0 


0.79 


1.4 


10.9 


0.0043 


75.7 


18.9 






(5.0 


1.05 


1.4 


14.2 


0.0054 


64.8 


16.2 






7.0 


1.46 


1.4 


16.7 


0.0071 


54.2 


13.6 


R2i = 


Phl 


4.0 


0.56 


1.4 


41.4 


0.038 


17.8 


4.5 






5.0 


0.79 


1.4 


71.0 


0.071 


11.7 


2.9 






(5.0 


1.05 


1.4 


109. 


0.12 


8.5 


2.1 






7.0 


1.46 


1.4 


153. 


0.20 


5.9 


1.5 



The initial system in each case consists of a neutron star (mass Mi) a donor at the onset of massive Case B evolution (mass M21, core 
mass Mzf). In the upper half of the table mass transfer is assumed to start when the donor reaches the Schonberg-Chandrasckhar limit, 
corresponding to the minimum possible orbital separation. In the lower half of the table mass transfer is assumed to start only when the 
donor has reached the Hayashi line, corresponding to the maximum possible orbital separation. The parameter (3qe = A Qce- 



tent rather than a transient LMXB (cf. King, Kolb & Bur- 
den 1996), so the system will eventually become transient. 
Once the envelope has been transferred, mass transfer will 
stop, and the system will become a pulsar binary with about 
the current orbital period P — 10 d, and a white dwarf com- 
panion with a mass which is slightly higher than that of the 
companion's present helium core. Clearly since the present 
core mass is at least 0.35Mq this P — combination will 
not obey the Hayashi-line relation (^]), so Cygnus X-2 will 
become another 'discrepant' system like those in Table 3. 

The reasoning of the last paragraph shows that Cygnus 
X-2 will cease to be a persistent X-ray binary within the 
current mass transfer time-scale tu = (Mi — M c )/M acc ~ 
10 7 yr. Its past lifetime as a persistent source before the cur- 
rent epoch, and its future one as a detectable transient after 
it, are both likely to be of a similar order, although full evo- 
lutionary calculations are required to check this. The fact 
that we nevertheless observe even one system like Cygnus 
X-2 strongly suggests that the birthrate of such systems 
must be relatively high, i.e. ~ 10~ 7 yr -1 in the Galaxy. 
Since the binary pulsar end-products of these systems have 
enormously long lifetimes, this may suggest that systems 
like Cygnus X-2 play a very important role in providing the 
Galactic population of millisecond pulsars. 

Cygnus X-2 thus fits naturally into a unified descrip- 
tion of long-period LMXBs in which super-Eddington Case 
B mass transfer is efficiently ejected by the neutron star. 
While the ejection process can already be inferred for the for- 
mation history of Hayashi-line LMXBs resulting from low- 
mass Case B evolution (Bhattacharya & van den Heuvel, 
1991; Kalogera & Webbink, 1996), Cygnus X-2 supplies the 
most powerful evidence that this process must occur. The 
work of Section 2 shows that it is very hard otherwise to 
reconcile the rather low current mass (M2 — 0.5 — 0.7Mq) 
of the secondary with its large radius (R2 — 7R©) and high 
luminosity (L2 — I5OL0). From Section 4 we see that the 
orbital period P — 9.84 d is far too long for the system to 
be the product of common-envelope evolution, leaving no 
realistic alternative for driving the required mass ejection. 
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8 APPENDIX: X-RAY HEATING IN CYGNUS 
X-2 

In Section 2.2 we considered a stripped subgiant model for 
Cyg X-2, and asserted that the observed orbital modulation 
of the optical flux (AV h B — 0.3 mag) would require an ex- 
tremely low inclination if one appeals to X-ray heating of 
the companion to raise its observed effective temperature to 
7330 K. Here we justify this claim. 

We consider a simple picture in which the hemisphere 
of the (spherical) companion facing the neutron star has ef- 
fective temperature 7330 K, while the other hemisphere has 
the Hayashi-line effective temperature 4100 K. We consider 
the effect of relaxing these assumptions below. Then viewing 
the heated face at the most favourable phase the observer 
sees hot and cool areas 27ri?|(l/2 + i/w), 2tt.R|(1/2 — i/w), 
where i is the inclination in radians, with the two expressions 
reversing at the least favourable phase. Neglecting limb- 
darkening, the ratio of maximum to minimum flux is 

F m ,x = (1/2 + yV)Bhot + (1/2 - i/7v)B cool 

F mi „ (1/2 + i/n)B cool + (1/2 - i/n)B hot ' 1 ' 

where Shot , Bcooi are the optical surface brightnesses of the 
hot and cool regions respectively. Approximating these by 
Planck functions at 5500 A, we find Bhot/S coo i ~ 15. Re- 
quiring Fmax/Fmin ^ 1.3 (AU b s 0.3 mag) in (|) shows 
that i 0.0745-7T, or i Ss 13.4°, as used in Section 2.2 

In reality the heated region would be smaller than a 
hemisphere, and its temperature higher than 7330 K in order 
to produce an average observed temperature of this value. 
However relaxing these limits clearly requires even smaller 
inclinations than the estimate above, because the contrast in 
optical surface brightness between the hot and cool regions 
would be even larger than the ratio ~ 15 we found above. 
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